1 



Lattice study of "fo(600) or a r 
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We investigate the propagator of "fo(600) or the a" by the full-QCD simulation with Wilson fermions. We 
calculate the mesonic correlator in the I — 0, J p = + channel on the 8 3 x 16 lattice. Plaquet action and Wilson 
fermion action are adopted. A coupling constant /3 is set to 4.8 and three kinds of hopping parameter, k=0.1846, 
0.1874 and 0.1891 are assayed. The disconnected diagram in the propagator is evaluated through taking average 
over 500 or 1000 Z2 noise. Simulations with the larger hopping parameter provide us with less noisy results. 
Though the statistics is not yet enough, our results indicate the existence of a pole with a mass in almost the 
same order as that of the p. 



1. Introduction 

In many text books, we can find the intuitive 
explanation of the spontaneous breaking of the 
the chiral symmetry based on the titled "Mexi- 
can Hat" potential, of which the angular and ra- 
dial modes correspond to the n and a mesons, re- 
spectively. Also in the the Nambu Jona-Lasinio 
model, the tr-degree of freedom plays an impor- 
tant role for the chiral symmetry breaking of 
QCD in the low-energy region. Despite its impor- 
tance, the experimental identification of the light 
scalar meson is a long standing puzzle and the 
a had disappeared from the list of Particle Data 
Group (PDG) for over 20 years. But, the scenes 
around the a are changing recently; re-analyses 
of 7T-7T scattering phase shift strongly suggest the 
existence of the I — and J PC = ++ meson 
and significant contributions of the a in the de- 
cay channels of heavy particles such as D— > tttttt 
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are reported [y|2| . The a may reveal itself also in 
T(35) -> Ttttt channel[||. In 1996 PDG ' "/o(400- 
1200) or c" appeared bellow 1 GeV mass region. 
In the 2002 edition, it appears as "/o(600) or cr" 
and as the 105th branch of the D + decay, "D 
— > a + 7r " is cited Q . Now it is an important 
task of the lattice calculation to confirm the prop- 
erties of the <7 meson based on QCD. 

In quenched QCD, DeTar and Kogut first mea- 
sured the a meson screening mass in a lattice sim- 
ulation || , Alford and Jaffe discussed the possible 
light scalar mesons as q 2 q 2 states ||, and masses 
and mixing of qq states and a glueball have been 
investigated by Lee and Weingarten 0. How- 
ever, as we have already reported in the previous 
proceedings ||, the connected and disconnected 
diagrams give almost the same amount of the con- 
tributions, indicating that the dynamical quark 
effect seems to be essential in the calculation. 

McNeile and Michael computed the mixed iso- 
singlet scalar masses of qq and glueball states in 
two kinds of situation, i.e, with and without the 
dynamical quark effects ||. We focus our atten- 
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tion to explore whether the a meson exists below 
1 GeV in this work. Based on the quite simi- 
lar motivation, Prelovsek presented in the present 
conference the calculation of the quenched scalar 
meson correlator with the domain wall fermions, 
where the disconnected diagram is taken into ac- 
count effectively through an approximate formula 
based on the chiral perturbation theory jl(J . 

2. Formulation of the a propagator 

As for the a meson, we adopt the operator as 



crlx 



y' y' M aOXO) + d c a (x)d c a (x) ^ 



C— 1 a— 1 



V2 



where u and d are the u-quark and d-quark Dirac 
spinors, respectively. The indices c and a de- 
note color and Dirac spinor indices, respectively. 
Quantum numbers of the operator are 1 = and 
J p = + . Because our simulation is full QCD, 
any kind of the state which has the same quan- 
tum number can contribute to the correlation as 
a virtual state. The a meson propagator is given 
by 

G{y,x) = (a(y)a(x)^) (y > x ) 

1 



Z 
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^2 u UyH(y) + d m d p(y) 

a, 6=1 a, 0=1 
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u%(x)u%(x) + d a a (x)d a a {x) 
V2 



t 



-S„-uWu-dWd 



= -(TrW~ 1 (x,y)W- 1 (y,x)) 

+2(TrW- 1 (y,y)TrW~\x,x)) . 

In Eq.(3), "Tr" represents summation over color 
and Dirac spinor indices. W^ _1 's are u and d 
quark propagators, U is the link variable of gluon, 
and S g is the pure gauge action of gluons. We 
assume that the u and d quark propagators are 
equivalent because u and d quark masses are al- 
most the same. From Eq.(3), we can see that a 
propagator is composed of the connected diagram 
(the first term of Eq.(3)) and the disconnected 
one (the second term of Eq.(3)). The quantum 
number of the a meson (1 = and J p = + ) is 
the same as that of the vacuum, and the vacuum 



(2) 
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expectation value of the a operator, (a(x)), does 
not vanish. Therefore, the contribution of (<j{x)) 
should be subtracted from the a operator. Thus, 



G(y,x) 



-(TrW- 1 (x,y)W- 1 (y,x)) 

-2(TrW / - 1 (y,y)TrW / - 1 (x,a;)) 

-2(Tr^- 1 ( 2/ ,y))(Tr^- 1 ( a :,x))(4) 



is adopted for the a propagator. The third term 
in Eq.(4) corresponds to the subtraction of the 
vacuum contribution. Our results show that the 
values of the second and the third terms in Eq. (4) 
are in the same order. Therefore, in order to ob- 
tain the signal correctly as the difference between 
these terms, the high precision numerical simula- 
tions and careful analyses are required. 

3. Results of numerical simulations 

We calculate the a propagator in the full QCD 
by using the Hybrid Monte Carlo (HMC) algo- 
rithm, in order to take the dynamical quarks 
which are important to estimate the disconnected 
diagrams of the a propagator. We use the Z 2 
noise method to calculate the disconnected di- 
agrams (the second and third terms in Eq.(|4|)). 
The 500 or 1000 random Zi numbers are gener- 
ated. The two-flavor Wilson fcrmion system is 
simulated on the 8 3 x 16 lattice. 

After 500 trajectories are updated in the 
quenched QCD, we start to update the configu- 
ration in the full QCD by using HMC algorithm. 
More than 500 trajectories by HMC are spent for 
thcrmalization. The a propagators are calculated 
on a configuration in every 10 trajectories. 

Based on ref.Q, we set (3 = 4.8 (a = 0.197(2) 
fm, k c = 0.19286(14)). Figures 1, 2 and 3 cor- 
respond to the three values of the hopping pa- 
rameter, k = 0.1846(fig. 1), 0.1874 (fig. 2) and 
0.1891(fig. 3), respectively. The numerical results 
are summarized in Table 1. 

4. Concluding Remarks 

We have reported our preliminary results on 
the a meson propagator based on a full QCD lat- 
tice calculation with dynamical fermions. The 
simulations with larger hopping parameters give 
us clearer signals. Though our simulation is still 
far from the chiral limit and no decay channel is 
yet open, our preliminary results obtained with 
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Table 1 

Summary of the results 




Figure 1. Propagators of ir, p and a for k 
0.1846 with periodic boundary condition. 




Figure 2. Propagators of ir, p and cr for n 
0.1874 with periodic boundary condition. 
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Figure 3. Propagators of ir, p and a for n 
0.1891 with periodic boundary condition. 



'),... , 

3 •'our result 



still a relatively low statics indicate the existence 
of the light a meson with the mass in almost the 
same order as that of the p meson. See Table 1. 
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